Orbital entanglement production in Andreev billiards with time-reversal symmetry 
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We study orbital entanglement production in a chaotic cavity connected to four single-channel 
normal-metal leads and one superconducting lead, assuming the presence of time-reversal symmetry 
(TRS). The scattered state of two incident electrons is written as the superposition of several two- 
outgoing quasi-particle states, four of which are orbitally entangled in a left-right bipartition. We 
calculate numerically the mean value of the squared norm of each scattered state's component, 
as functions of the number of channels in the superconducting lead. Its behavior is explained as 
resulting from the proximity effect. We also study statistically the amount of entanglement carried 
by each pair of outgoing quasi-particles. When the influence of the superconductor is more intense, 
the average entanglement is found to be considerably larger than that obtained using normal cavities. 

PACS numbers: 73.23.-b,03.67.Bg,05.45.Mt 



I. INTRODUCTION 



Quantum correlations are in the very heart of quan- 
tum physics^ and of applications like quantum computing 
and quantum cryptography^ The potential of quantum- 
transport devices to produce and manipulate entangle- 
ment has been explored during the last years.^ Several 
of those devices produce entanglement from scattering 
processes taking place in its components. 

Chaotic cavities have been proposed as orbital 
cntanglers£ The cavity is connected to four one-channel 
normal-metal leads, two at the left and two at the right. 
An electron leaving the cavity at either the left or the 
right represents a qubit, and entanglement between two 
qubits can be studied considering a left-right biparti- 
tion. The mean value and the variance of the concur- 
rence, a quantifier of entanglement, were initially calcu- 
lated in the context of random matrix theory^ followed 
by more complete statistical analysis.— Non-ideal con- 
tacts have also been considered^ Constraints on en- 
tanglement production imposed by the geometry of the 
device were explored^ and it was found that more en- 
tangled states are less likely to be produced in general. 

Scattering processes in normal-superconducting (NS) 
hybrid systems are another source of entanglement ^LiS 
Entangled pairs of quasi-particles can be generated af- 
ter Andreev reflections which take place at the normal- 
superconductor interface. In the present work, we pro- 
pose the use of Andreev billiards as orbital entanglers. 
We study the generation probability of each pair of quasi- 
particles, and how this quantity and the amount of en- 
tanglement are affected by the influence of the super- 
conductor. We find a notable increase in the production 
of entanglement induced by the proximity effect ^1 when 
compared to normal cavities. On the other hand, we ob- 
serve in general a compensation between the amount of 
entanglement and the probability of producing it. 

The paper is organized as follows. In Sec. |TT]we explain 
how the device is designed, and what is the structure of 
its scattering matrix. We also write the scattered state as 
a function of the transmission properties of the system. 
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FIG. 1: (Color online) A chaotic cavity is attached to four 
normal leads, denoted by L\, Li, R\ and R2, and one floating 
superconductor lead. Each normal lead has one conduction 
channel while there are N 3 open channels in the superconduc- 
tor lead. A potential V is applied on the left leads and the 
right leads are grounded. 



The mean value of the squared norms of states describing 
two-outgoing quasi-particles, even those non-entangled, 
is analyzed in Sec. IIIII Statistics of concurrence is stud- 
ied in Sec. HVl We present its full distribution, as well 
as its mean value and its variance, for different values of 
the number of open channels in the superconducting lead. 
We summarize and conclude in Sec. [Vj Two additional 
sections are included as appendices. The equations used 
to perform the numerical simulations arc presented in 
Appendix [Al while some analytical considerations about 
the squared norms of the entangled states appear in Ap- 
pendix [B] 



II. DESIGN OF THE DEVICE 

The setup is represented in Fig. Q] A central chaotic 
cavity is connected to five leads. Four leads are normal 
and have only one open conduction channel; they are 
represented at the left and right sides of the cavity and 



denoted by L±, L2, Ri, and i?2- A superconducting lead 
with N s channels is attached to the top of the cavity. 
A potential V is applied on the left leads and the right 
leads are grounded. The superconducting lead is float- 
ing. Thus an electrical current crosses the device by the 
normal leads. Specifically, pairs of electrons with ener- 
gies between the Fermi energy Ep and Ep + eV enter the 
cavity from the left leads. It is assumed that the energy 
of the incident electrons is much smaller than the gap A 
in the superconductor, i.e. eV <C A. This implies that 
charge transfer in the NS interface takes place only via 
Andreev reflection^ Since the device has only one su- 
perconducting lead, the phase of the superconductor is 
irrelevant. We also neglect temperature fluctuations. 



A. Scattering matrix of the system 

The potential V is small enough to neglect the depen- 
dence on energy of the scattering matrix. We emphasize 
that this matrix does not belong to any universal class 
of random matrices. In order to reach universality in 
the presence of TRS, it would be necessary to attach an- 
other superconducting lead to the cavity, with a phase 
difference of ttJ^^- Also the number of open channels 
in the superconducting lead must be larger than oneJ£ 
We focus on a different regime, where proximity effects 
are important. The scattering matrix will be denoted by 
Sns and has the following structure: 



Sns 



f t' 

i f 



(i) 



The blocks f if') and t (f) contain, respectively, the 
reflection and the transmission amplitudes for quasi- 
particles coming from the left (right). Each block has 
its own electron-hole structure, which is indicated with 
the hat symbol. For example, the matrix f is 



(2) 



On the other hand, each sub-block r a $ in the last equa- 
tion is a 2 x 2 matrix, whose elements are the reflection 
amplitudes for the leads L\ and L^. They are given by 



r.a/3 



a/3 aB 

r ll r 12 

a0 a(S 

r 21 r 22 



(3) 



The other blocks of Sns have the same structure. The 
scattering matrix of the whole system can be expressed 
as a function of the scattering matrix of the normal 
cavity^ which belongs to the Circular Orthogonal En- 
semble. Equations showing this dependence are pre- 
sented in Appendix lAl 



B. Scattered state 

The scattered state of two incident electrons can be ex- 
pressed as |* SCBt ) = |*LL> + |*fii?.) + |*Lfl>- State |* LL ) 
(\^rr)) characterizes two quasi-particles being scattered 
to the left (right), so is not entangled. One quasi-particle 
scattered to the left and the other one to the right is rep- 
resented by the state \^lr)- In order to give explicit 
expressions for these three states, we define the projec- 
tor matrix 



1 




(4) 



The structure of this matrix accounts for the fact that 
the incident state is given by two electrons coming from 
the left. It is also convenient to define the operators £" 
and 1Z" , which create a quasi-particle of type a going 
out at the lead Lj and Rj, respectively, with energy e if 
a = e and — e if a = h (Ep < e < Ep + eV). We also 
define 



5-- 

V 
7 : 



ra r , 
= ta'i T , 
= ra't T , 



(5) 
(6) 

(7) 



where the superscript "T" indicates the transposition op- 
eration. Using the above objects, the three states can be 
expressed in the following way: 



i*">4E^ Q£ "i°>< 



jk-a/3 



i*«fl> = 5E^^i^i°)' 



^lr) = 



jk;a/z 

jk;a.f3 



7^7^|0>, 



(8) 

(9) 
(10) 



where |0) is the state without electronic excitations. Each 
term in these sums represents a pair of quasi-particles of 
types a and (3 getting out at the leads j and k. The 
equations presented in this section are general and do 
not specify any particular property of the NS system. 



III. STATISTICS OF THE SQUARED NORMS 

In this section wc focus on the probability of producing, 
during a fixed time, a given number of pairs a and (3 being 
scattered to the left, to the right, or to different sides. 
This probability is determined by the squared norm of the 
state which represents each process. Let us denote such 



squared norms as A/"™ , N^ , and J\f a " . We are assuming 
that A/" a/3 refers to the quasi-particle a scattered to the 
left and the other one, /3, to the right. For the other 
two squared norms there is no distinction between a, f3. 
Using Eqs. (|51)- (fTU)l . they are expressed as functions of 
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FIG. 2: (Color online) Mean values of the squared norms 
A/" , Af£ , and Af^ are represented in panels a), b), and 
c) respectively, for different numbers of open channels in the 
superconducting lead. 10 5 samples were generated for each 
value of N a . 



the transmission properties of the system through the 
following equations: 



Aff = Tr(S a0 S a ^) 



Tr(77 Q/3 77 Q/3t ) 
Af afi = Tr( 7 Q ^7 Q/3t ) 



Aff 

r<xB 



(11) 
(12) 

(13) 



Let us discuss the underlying mechanism behind the 
influence of proximity effects on transport properties. As 
the number of open channels in the superconducting lead 
increases, the probability of incident electrons to undergo 
an Andrcev reflection in the NS interface and thus to 
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FIG. 3: Function E r (in percent) for values of iV s varying 
from 10 to 500. Each dot was calculated using 5 x 10 3 random 
matrices. 



be backscattered as holes to the left before "feeling" the 
chaos of the cavity gets larger. As N s — > oo this probabil- 
ity becomes unit, and the device works as if all electrons 
were Andreev reflected upon arrival at the cavityJ^ 

A consequence of direct Andreev backs cattering is an 
increased number of holes leaving the cavity at the left, 
and therefore a decreased transmission of any quasi- 
particle to the right. This is verified in Fig. [2j As N s 

gets bigger, AT , AT , J7 L , AT fl , and J7 R monotonically 
decrease, because the probability of an electron to be 
scattered to the left continually decreases. On the other 

hand, AT , AT" and A~f L have a non-monotonic behav- 
ior, initially increasing due to direct Andreev backscat- 

tcring. The quantity Af decreases more slowly than the 
other squared norms of entangled states because scatter- 
ing to the right is more probable for electrons than for 

holes. Furthermore, AT L (A^_r ) monotonically increases 
(decreases), also because of direct Andreev backscatter- 
ing. 

The dominance of A/" over the other entangled com- 
ponents can be quantified by 



E r = 



AT+ A7 e 



W 



hi, 



Af 



he 



(14) 



The behavior of E r is shown in Fig. [31 where N s runs 
from 10 to 500. Notice that for N s = 500, the sum of 

the other components represents less than 5% of Af . 
Although all four kinds of events are rare for large values 
of N Sl a setup for detection of entanglement would char- 
acterize correlations between one left-outgoing hole and 
one right-outgoing electron with a small relative error. 

Direct Andreev backscattering docs not occur if TRS 
is broken, because a random phase is accumulated dur- 
ing the process and so it is killed after averaging. An 
implication of this is that in the presence of a magnetic 
field and for N s S> 1, the statistical distributions of Af a ^ 
are very similar for any a and /3. We corroborated this 
observation through numerical simulations (not shown) . 
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FIG. 4: (Color online) Distributions of concurrence are rep- 
resented for N s equal to a) 5, b) 25, and c) 100. Each distri- 
bution was made using 10 6 random matrices. 



IV. STATISTICS OF CONCURRENCE 



The concurrence of the entangled states is given b y 19 i 20 



C ali = 2 V /Det( 7 a 'V/ 3t )/Tr(7 Q 'V' 3t ). (15) 

In contrast with the result obtained for a generic nor- 
mal system^ this expression depends not only on the 
eigenvalues of the normal scattering matrix, but on its 
eigenvectors as well. This prevents us from proving the 
geometric constraint for the sq uared norm and the con- 
currence discovered in Ref. [[(J for normal systems, that 
forces 7V(1 + C) < 1. Nevertheless, we have observed in 
numerical simulations that this constraint seems to re- 
main valid for the considered system. Using Eqs. (J7J, 
(I15p . and taking into account the unitarity of Sns, a 
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FIG. 5: (Color online) a) Mean value and b) variance of the 
concurrence calculated from 10 random matrices for each 
value of N s . 



weaker constraint can be proved, 

N al3 C a(5 < 1/2. (16) 

This inequality is not related to specific physical prop- 
erties of the system, but it is an implication of unitarity 
and the design of the device. 

Fig. 2] shows the distributions of concurrence for 
N s = 5, 25, and 100. They suggest that all pairs of left- 
right outgoing quasi-particlcs are entangled, but none of 
them is maximally entangled. It is clearly seen that C he 
is the quantity most sensitive to the presence of the su- 
perconductor. However, variations are small for all the 
distributions when N s is large, which indicates conver- 
gence, as can be observed in Fig. [5] where we show the 
average value and the variance of the concurrence. 

As the number N s increases, there is a remarkable 
enhancement in the production of entanglement for the 

—he 

hole-electron component: C reaches values above 0.58, 
reflecting a considerable increase compared to the nor- 
mal case, where C rj 0.38£ Note that the other three 
components have only slight dependence on N s . More- 
over, variations of C he from sample to sample are more 
concentrated around its mean value than for the other 
components. The variance of this quantity decreases as 
N s increases (see Fig. \5}p). We found that var(C ) os- 
cillates around 0.53 for large N s . This value is slightly 
less than the variance for the other components, as well 
as those found for a normal cavity^ 



SUMMARY AND CONCLUSIONS 



and 



We studied orbital entanglement production in An- 
dreev billiards, assuming time-reversal symmetry in the 
system. Numerical simulations were performed to statis- 
tically analyze the squared norm and the concurrence of 
states describing a pair of outgoing quasi-particles. The 
behavior of the squared norm for varying numbers of 
channels in the superconducting lead was explained as 
a consequence of Andreev reflection. Even though the 
norms of the entangled components decrease as the in- 
fluence of the superconductor is stronger, the squared 
norm of the state describing one left-outgoing hole and 
one right-outgoing electron always dominates. In this 
regime, there is a notable enhance in the amount of en- 
tanglement for that state. The mean concurrence reaches 
values greater than 0.58, exceeding the value found for a 
normal cavity. 

There exist a general feature in the entanglement pro- 
duction using the geometry studied in this paper and 
other works£~— The proht obtained in the amount of en- 
tanglement is, in some way, compensated by a decreased 
probability of producing the state. In the case of the sys- 
tem under consideration, this tendency is reinforced by 
the proximity effect. 

It is still an open problem to establish a suitable setup 
for detection of the entanglement carried by each type of 
outgoing pair of quasi-particles. 







Vri Ns . 



(A4) 



r being the transparency of the barrier, and I n the nx n 
identity matrix. 

First, it is necessary to compose So and S c , to obtain 
the total normal scattering matrix Sn ■ This is done using 
the following equations: 



rjv = r c + t c M N r Q t c , 
tN = i to (in + r c M N r )t c , 
r' N = -r' - t r c M N t' , 
t N = i t c MNt , 



(A5) 

(A6) 
(A7) 
(A8) 



where M^ = (Jjv s — ^o r c) _1 - The phase of the super- 
conductor is irrelevant and it can be set equal to zero. 
Neglecting the energy-dependence, the blocks of Sns in 
the electron-hole space can be written as^ 



S% s = r N -t' N (r' N )*M e t N , (A9) 

-i (t' N )*M e t N , (A10) 



Qeh 



where M e = (I Ns +r' N (r' N )*) 
is structured as 



On the other hand, S^f 
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S^g and S$g are not necessary to calculate norm and 
concurrence, because we consider the initial state given 
by two electrons entering the cavity from the left. 



Appendix A: Equations to perform the numerical 
simulations 



Appendix B: Some analytical considerations about 



Let's define the matrix of the normal cavity as 

s °=0o°. ! !)' <A1) 

Our numerical simulations arc performed through ran- 
dom generations of So, which has dimensions 2N S x 2N S . 
We need to close N s — 4 conduction channels in the nor- 
mal side (assuming N s > 4), which is done introducing a 
fictitious barrier with very small transparency. Only four 
channels remain open, one for each normal lead. The 
scattering matrix for the fictitious connector is given by 



S c 



t' 



We choose r c = r',t c = t', and define 



r..r^ 



V x/T - Tl 



iVa-4 



(A2) 



(A3) 



We proceed analogously to Ref. |2lj, where a similar 
device was studied. However, an essential difference is 
the small numbers of open channels in the normal leads, 
which requires to take into account all diagrams. 

It is convenient to parameterize the scattering matrix 
of the cavity through the polar decompositio n 18 ! 22 



So 



/ Vyfh^Tv 1 






v(0 iy/T) 



\ 







Vh^T 



) 

(Bl) 
In the above equation v and w have dimensions 4x4 and 
N s x N s , respectively, and belong to unitary circular en- 
sembles. They account for the isotropy assumption.— On 
the other hand, T = diag(T^ 1 , tx 2 , trj , tr 2 ) is a function 
of the transmission eigenvalues. The eigenvalue T£.(_r.) 
quantifies the probability that an electron entering from 
lead Lj (Rj) meets the superconductor before leaving the 
cavity. 
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FIG. 6: The squared norms of the entangled components are 
sketched using the symbology of the diagrammatical method 
for averaging over the unitary group. 



Making S N = So in Eqs. ([Ml) and (|A10j) . and using 
the polar decomposition, the equations for the squared 
norms of the entangled components can be written as 

N ee = -Tr [P L vT ee v T avT ee v T P R v*T ee ^ av*T ee ^} , 

(B2) 

N eh = -Tr [P L vT ee v 7 'avT^v 1 P R vT he v ] av*T ee v^] , 

(B3) 

N he = -Tr [P L v*T he v T <jvT ee v T P R v*T ee v^av*T he v T ] , 

(B4) 

N hh = -Tr [P L v*T he v T <7vT he v^P R vT he v f av*T he v T ] , 

(B5) 
where we have defined the projectors 

.(\ 
.0 



(J oj 



Expressions (|B2|) - (|B5[) can be represented, using the 
symbology of the diagrammatical method for averaging 
over the unitary group,— as shown in Fig. |6] Each part of 
the representation indicates some kind of event. Matrices 
T ae connected to two dashed lines refer to a scattering 
process in the NS system, meaning that an electron enters 
and a outs. On the other hand, Pl (P R ) selects processes 
describing a quasi-particlc going out at the left (right), 
and a assures that two electrons enter from the left. 

In order to make the averaging, all topologies result- 
ing from matching white/black dots should be generated. 
There arc 24 x 24 = 576 different topologies for each 
quantity. Every diagram contains 7 matrices in its top, 
indicating a process in which two electrons enter the cav- 
ity from the left (matrix a), one is scattered to the left 
(matrices are connected to Pl) and the other one to the 
right (matrices are connected to P R ) . The bottom repre- 
sents the same process, but inverted in time. Note that 
the chronological order of the events is determined by the 
sequence of matrices, and not by the arrows. 

T-cycles generated from connecting black dots indicate 
that most of the diagrams are zero. If the black dots 
annexed to P R are matched to any other black dot, a 
null diagram is generated. This happens because two 
events taking place at different sides of the device can 
not be directly connected. For this reason, the number 
of diagrams different from zero is 144, which is still large. 

Even though it is hard to calculate all the diagrams, it 
is possible to know their dependence on the transmission 
eigenvalues, which results from T-cycles generated from 
white dots. There is indeed a small number of possible 
functional dependences for each squared norm, as shown 
bellow 



N e 



(1- 7-0(1 -Tj) 

(2-r 4 ) 2 (2-r,)2 



A' 



hh 



feh . 


(1 - Ti)TjT k 




1 ' (2- 


-r l ) 2 (2-r j )(2- 

TiTjTkTl 


- Tk)' 



(2-r0 2 (2-r,)(2-r fe )(2-r0' 



(B9) 



(BIO) 



(Bll) 



We have assumed that 1 < i,j, k, I < 4. The dependence 
for J\f he is the same as the one for J\f eh . 

The statistical distribution of the transmission eigen- 
values is given b y 24 i 25 



T ee and T he , on the other hand, encode the transmission 
properties of the NS system. They are given by 



T ee = (2y/l 4 -T)(2I 4 - T)-\ (B7) 

T he = -T(2h-T)- 1 . (B8) 



P(r)=C [J iTi-TiKwmTjW'-W, (B12) 

l<i<j<4 

where C is a normalization constant. For increasing N s , 
this distribution is negligible except around the point 
n = 1. Considering the scaling T{ = l — £i/N s , it is possi- 
ble to find that for large N s the mean value J\f ~ l/-^ 2 - 
The other components are proportional to 1/N S . In this 



analysis, we have neglected terms of higher order in N s 1 , 

and assumed that AT is zero in the limit N s — > oo. 

It is also possible to express the squared norms as a 
function of the charge cumulants of a normal device. As 
an example, let's do this for J\f ee . Making r^ + l/r^ ps 2 
and using the Levitov-Lesovic formula^ we find that 



[TT(T ee ) 2 } 2 ^^2-q3) 2 , 



(B13) 



We have denoted with qj the cumulant of transmitted 
charge of order j, for a normal device consisting of a 
chaotic cavity connected to two leads with 4 and N s 
open channels. Following the same argument, concur- 
rence might also be written as a function of the charge 
cumulants. 



and 



Tr(T e 



2[6g 2 - 5(<7 3 + Q4- q 5 ) - <fe]/15. (B14) 
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